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ABSTRACT 

We investigate a stationary pair production cascade in the outer magnetosphere of 
^ ' a spinning neutron star. The charge depletion due to global flows of charged particles, 

p causes a large electric fleld along the magnetic fleld lines. Migratory electrons and/or 

Q ! positrons are accelerated by this fleld to radiate curvature gamma-rays, some of which 

^ I collide with the X-rays to materialize as pairs in the gap. The replenished charges 

^ ■ partially screen the electric fleld, which is self-consistently solved, together with the 

distribution functions of particles and gamma-rays. By solving the Vlasov equations 
describing this pair production cascade, we demonstrate the existence of a stationary 
^ I gap in the outer magnetosphere of PSR B1055-52 for a wide range of current densities 

■ flowing in the accelerator: From sub to super Goldreich- Julian values. However, we 

flnd that the expected GeV spectrum becomes very soft if the current density exceeds 
the Goldreich- Julian value. We also demonstrate that the GeV spectrum softens with 
decreasing magnetic inclination and with increasing distance to this pulsar. We thus 
conclude that a sub-Goldreich- Julian current, a large magnetic inclination, and a small 
distance (500 pc, say) are plausible to account for EGRET observations. Furthermore, 
it is found that the TeV flux due to inverse Compton scatterings of infrared photons 
whose spectrum is inferred from the Ray leigh- Jeans side of the soft blackbody spectrum 
is much less than the observational upper limit. 



Subject headings: gamma-rays: observations - gamma-rays: theory - magnetic flelds - 
pulsars: individual (B 1055-52) 
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1. Introduction 

The EGRET experiment on the Compton Gamma 
Ray Observatory has detected pulsed signals from 
seven rotation-powered pulsars (for Crab, Nolan et 
al. 1993, Pierre et al. 1998; for Vela, Kanbach et 
al. 1994, Fierro et al. 1998; for Geminga, Mayer- 
Hasselwander et al. 1994, Fierro et al. 1998; for PSR 
B1706-44, Thompson et al. 1996; for PSR B1046-58, 
Kaspi at al. 2000; for PSR B1055 52, Thompson et 
al. 1999; for PSR B1951+32, Ramanamurthy et al. 
1995). The modulation of the 7-ray hght curves at 
GeV energies testifies to the production of 7-ray radi- 
ation in the pulsar magnetospheres either at the polar 
cap (Harding, Tademaru, & Esposito 1978; Daugh- 
erty & Harding 1982, 1996; Michel 1991; Sturner, 
Dermer, & Michel 1995; Shibata, Miyazaki, & Taka- 
hara 1998; Harding & Muslimov 1998; Zhang & Hard- 
ing 2000), or at the vacuum gaps in the outer mag- 
netosphere (e.g., Cheng, Ho, & Ruderman 1986a,b, 
hereafter CHRa,b). Both of these pictures have had 
some success in reproducing global properties of the 
observed emission. However, there is an important 
difference between these two models: An polar gap ac- 
celerator releases very little angular momenta, while 
outer gap one could radiate them efficiently. More 
specifically, the total angular momentum loss rate 
must equal the energy loss rate divided by the angu- 
lar velocity of the star, implying an average location 
of energy loss on the light cylinder (Cohen & Treves 
1972; Holloway 1977; Shibata 1995), of which distance 
from the rotation axis is given by 

t^LC = ^, (1) 

where il denotes the angular frequency of the neutron 
star, and c the speed of light. 

On these grounds, the purpose here is to explore 
a little^ further into a model of an outer gap, which 
is located near the light cylinder. If the outer mag- 
netosphere is filled with a plasma so that the space 
charge density is equal to the Goldreich-Julian den- 
sity, pgj = — fi-Bz/(27rc), then the field-aligned elec- 
tric field vanishes, where is the component of the 
magnetic field along the rotational axis. However, the 
depletion of charge in the Goldreich and Julian model 
in a region where it could not be resupplied, may 
cause a vacuum region to develop. Holloway (1973) 
pointed out the possibility that a region which lacks 
plasma is formed around the surface on which pcj 
changes its sign. Subsequently, CHRa,b developed 



a version of an outer magnetospheric 7-ray emission 
zone in which acceleration in the Holloway gaps above 
the null surface, where vanishes, brought particles 
to large Lorentz factors (~ 10^'^). 

After CHRa,b, many outer-gap models have been 
proposed (Chiang & Romani 1992, 1994; Romani 
& Yadigaroglu 1995; Romani 1996; Zhang & Cheng 
1997; Cheng, Ruderman & Zhang 2000). Their basic 
idea is as follows: A deviation of the charge density 
from pGJ results in an electric field along the magnetic 
field, B. If this electric field becomes strong enough to 
accelerate pairs to ultrarelativistic energies, they 
could radiate 7-rays tangential to the curved B field 
lines there. The curvature 7-rays may be converted 
into pairs via two photon collisions. In order to 
keep a steady current flow and the charge density poj 
in the regions outside the gap, the gap will grow un- 
til it is large enough and the electric field is strong 
enough to maintain a sufficient supply of charges to 
the rest of the open field line region. If the gap ends 
in a region where pG,i 7^ and E - B 7^ 0, charges from 
the surrounding region may flow in through the end 
(or the boundary). If both boundaries are located on 
the null surface, pairs produced in the gap will re- 
place the charge deficiency inside the gap, and finally 
the gap will be filled up. However, if a vacuum gap 
extends to the light cylinder, the charged particles 
created in the gap should escape from the magneto- 
sphere, so the gap would not be quenched. Hence, 
stable outer gaps (if they exist) are those from the 
null surface to the light cylinder along the last-open 
field lines. In an outer gap, the inner boundary lies 
near the intersection of the null surface where pcj 
vanishes and the boundary of the closed field lines of 
the star on which the magnetospheric current does 
not exits. 

Recently, solving the Vlasov equations that de- 
scribe a stationary pair production cascade, Hirotani 
and Shibata (2001a,b, Papers VII and VIII) eluci- 
dated important characteristics of the particle ac- 
celeration zones in the obiter magnetosphere. They 
demonstrated that a stationary gap does not ex- 
tend between the null surface and the light cylin- 
der. Rather, their width (W) along the field lines 
is adjusted so that a single particle emits copious 7- 
rays one of which materialize as a pair in the gap 
on average. The resultant W becomes about 5% of 
n7LC in the case of the Crab pulsar. The produced 
pairs in the gap do not completely cancel the accel- 
eration field, because the particles do not accumulate 
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at the gap boundaries. Outside of the gap, the par- 
ticles flows away from the gap along the magnetic 
field lines as a part of the global current flow pat- 
tern. In another word, the remained, small-amplitude 
electric field along the magnetic field lines outside of 
the gap prevents the gap from quenching. Moreover, 
they demonstrated that the gap position is not fixed 
as considered in previous outer gap models. Their 
position shifts as the magnetospheric current flowing 
the gap changes. That is, if there is no particle injec- 
tion across either of the boundaries, the gap is located 
around the null surface, as demonstrated in their ear- 
lier papers (Hirotani & Shibata 1999a, b,c, Papers I, 
II, III; Hirotani 2000a, b,c, Papers IV,V,VI; see also 
Beskin et al. 1992 and Hirotani & Okamoto 1998 for 
a pair-production avalanche in a black-hole magne- 
tosphere). However, the gap position shifts towards 
the light cylinder (or the star surface) if the injected 
particle flux across the inner (or the outer) bound- 
ary approaches the Goldreich- Julian value. In other 
words, the accelerator can be located at any place in 
the magnetosphere; their position is primarily con- 
strained by the magnetospheric current. 

In Papers VII and VIII, they examined the outer 
gaps formed in the magnetospheres of young or mil- 
lisecond pulsars, assuming a sub-Goldrcich- Julian par- 
ticle injection across the boundaries. In this paper, 
we apply the same method to a middle-aged pulsar 
B1055-52 and demonstrate that a stationary gap is 
maintained even if the injected particle flux exceeds 
the Goldreich-Julian value, and that the possibility of 
such a large particle flux injection is ruled out if we 
compare the predicted 7-ray spectrum with observa- 
tions. 

In the next section, we formulate the basic equa- 
tions describing the stationary pair production cas- 
cade in the outer magnetosphere. In § 3, we ap- 
ply the method to a middle-aged pulsar, B 1055-52 
and reveal that large magnetic inclination and a sub- 
Goldreich- Julian current are plausible for the outer- 
magnetospheric emission from this pulsar. In the flnal 
section, we compare the results with previous works. 

2. Basic Equations 

We first consider the continuity equations of par- 
ticles in § 2.1, the Boltzmann equations of 7-rays in 
§ 2.2. Combining these equations with the Poisson 



accelerator along the field lines in § 2.3. We finally 
give appropriate boundary conditions in 



2.4 



equation, we present the Vlasov equations describing 
the one-dimensional electrodynamic structure of an 



2.1. Particle Continuity Equations 

In the same manner as in Papers I- VIII, we sim- 
ply assume that the electrostatic and the curvature- 
radiation-reaction forces cancel each other in the 
Boltzmann equations of particles. Under this mono- 
energetic approximation, the Boltzmann equations 
reduce to continuity equations. Assuming that the 
toroidal bending is negligible for the magnetic field, 
we obtain the following stationary {dt + = 0) 
continuity equations (Paper VII), 



±B- 



ds \ B 



1 

c Jo 



(2) 



where 7V_|_ and N_ represent the spatial densities 
of positrons and electrons, respectively, G+ and G_ 
the distribution functions of outwardly and inwardly 
propagating 7-rays, respectively, and B the strength 
of the magnetic field. The pair production redistribu- 
tion functions r]p± can be defined as 



Vp±i^'r) = (1 - Mc)c 



rfex^^crp(e^,ex,/ic), (3) 



where Cp is the pair-production cross section and 
cos~^/ic refers to the collision angle between the 7- 
rays and the X-rays with energy rricC^e^. X-ray num- 
ber density between dimensionless energies ex and 
Ex + de^, is integrated over Cx above the threshold 
energy eth = '2(1 ~ iJ,c)^^f-^^ to infinity. 

In Paper VII, we adopted — cos(T/F/ti7Lc) for 
the Crab pulsar magnetosphere, in which a power-law 
X-ray component dominates the surface ones. How- 
ever, in the present paper, we apply the theory to 
a middle-aged pulsar, B1055-52, of which soft X-ray 
emission is dominated by surface components. We 
thus evaluate fic by the collision angle between the 
surface X-rays and the curvature 7-rays. The colli- 
sion approaches head-on (or tail-on) for inwardly (or 
outwardly) propagating 7-rays, as the gap shifts to- 
wards the star. 

2.2. Boltzmann Equations for Gamma- rays 

For simplicity, we assume that the outwardly (or 
inwardly) propagating 7-rays dilate (or constrict) at 
the same rate with the magnetic field; this assumption 
gives a good estimate when W ^ zuj^c holds. Under 
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this assumption, we obtain (Paper VII) 
OS \ B c c 



(4) 



where the redistribution function, 77c, for the curvature- 
radiation is expHcitly defined in equation (12) in Pa- 
per VII (see also Rybicki, Lightman 1979). 

2.3. One-dimensional Analysis 

As described at the end of § 3 in Paper VII, it is 
convenient to introduce the Debey scale length c/wp, 
where 



' 47re2 VIB, 



27rce ' 



(5) 



and Bent is the magnetic field strength at the gap 
center. The dimensionless coordinate variable then 
becomes 

e ^ K/c)s. (6) 

By using such dimensionless quantities, and by as- 
suming that the transfield thickness of the gap is 
greater than W , we can rewrite the Poisson equation 
into the following one-dimensional form (Paper VII): 



-Ell = 



dE 



[n+iO-n-iOl + 



(7) 



(8) 



Bent Bent 

where = e'^{s)/ (mcC^); the particle densities are 
normalized by the Goldreich- Julian value as 



2'Kce N. 



± 



n B 



(9) 



We evaluate B^/B at each point along the last-open 
field line, by using the Newtonian dipole field. 

Let us introduce the following dimensionless 7-ray 
densities in the dimensionless energy interval between 
and Pf. 



2'Kce 



VLB, 



cnt J/3i_ 



de^G±{s, e^). 



(10) 



In this paper, we set /3o = 10^, which corresponds to 
the lowest 7-ray energy, 51.1 MeV. We divide the 7- 
ray spectra into 11 energy bins and put /3i = 10^ '^, 
132 = 103, /33 = 103■'^ Pi = 103■^ 05 = 10=^■'^ Pe = 



J2 = -LU , P3 
10* 1^- — in4-25 



Pr = 104■2^ Ps = lO^-s, Pg = 104-75^ p-_ 
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10^, Pii = lO^'^^. We can now rewrite the continuity 



quation 



of particles into 
9 

B ' 



^ = ±%E[WffV(?)+^p-'ff-(e)], (11) 



where the magnetic field strength, B, is evaluated at 
each ^. The dimensionless redistribution functions 
rjpj. are evaluated at the central energy in each bin as 



1 



-Vp± 



Pi-i+Pi 



(12) 



A combination of equations (|l l| ) gives the current 
conservation law. 



jtot = "-+(0 + "--(0 — constant for ^. 



(13) 



If itot = 1 holds, the current density per unit mag- 
netic flux tube equals the Goldreich- Julian value, 
17/(27r). 

The Boltzmann equations (^) for the 7-rays are in- 
tegrated over between dimensionless energies Pi-i 
and Pi to become 



d i d , i , I vlBjO 

-g^ = - (In B)g^TV,± 9^±^^n^, 



where z = 1, 2, • • •, m (m = 9) and 



(14) 



du / K^{t)dt (15) 



is dimensionless; ^^5/3, i?c and F refer to the modified 
Bessel function of 5/3 order, curvature radius of the 
magnetic field lines, and the particle Lorentz factor, 
respectively. 

To ensure that the Lorentz factors do not exceed 
the maximum attainable limit, ^'(^2)1 we compute P 
with (Paper VII) 



1 
f 



1 



1 



^2(6)' 



where 



V 2e 



d^ 



ds 



1 



1/4 



(16) 



(17) 



represents the terminal Lorentz factor that is realized 
if the electric force, e\d'^ /dx\, balances with the radi- 
ation reaction force. 
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2.4. Boundary Conditions 

We now consider the boundary conditions to solve 
the Vlasov equations (0), (|), (0), and (|l|). At the 
inner (starward) boundary {S, — 'we impose (Pa- 
per VII) 

E\\i^i) = 0, (18) 

v(a) = 0, (19) 

ffV(a)=0 (z-l,2,...,9), (20) 

n+{^,)=ji, (21) 

where the dimensionless positronic current at ^ = ^ 
is denoted as ji. Note that equations ( p^ and ( pl| ) 
yield 

^-(6)=Jtot-Ji. (22) 
At the outer boundary (^ = a), we impose 

^||(6) = 0, (23) 

ffL(6) = (z-l,2,...,9), (24) 
ri~(6)=j2. (25) 

The current density created in the gap per unit flux 
tube can be expressed as 



Jgap — Jtot ^ Ji ^ J2- 



(26) 



We adopt jgap, ji, and j2 as the free parameters. 

We have totally 2m + 6 boundary conditions (p^)- 
(25) for 2m + 4 unknown functions \E', n±, g\. 
[i = 1,2,- • -jto), where m = 11. Thus two extra 
boundary conditions must be compensated by making 
the positions of the boundaries ^ and ^ be free. The 
two free boundaries appear because E\\ = is imposed 
at hoth the boundaries and because jgap is externally 
imposed (see § 4T for details). In other words, the 
gap boundaries (^ and a) shift, if ji and/or varies. 

It is worth mentioning that the gap width, W , is 
related with the X-ray field density as follows (for 
details, see Papers V, VII, and VIII): 



W 



Ap jgap 
jtot ' 



(27) 



where Ap refers to the pair production optical depth, 
and the expectation value of the 7-ray photons 
emitted by a single particle while it runs the gap. 
Equation (|^ is automatically satisfied by the sta- 
tionary Vlasov equations. From equation (p7|), we 
can expect an extended gap for middle-aged pulsars, 



because their less dense X-ray field leads to a large Ap, 
and hence a large VF, compared with young pulsars 
like Crab. We examine such features more accurately 
in the next section, by analyzing the Vlasov equations 
numerically. 

3. Application to PSR B1055-52 
3.1. X-ray and Infrared Field 

Combining ROSAT and ASCA data, Greiveldinger 
et al. (1996) reported that the X-ray spectrum con- 
sists of two components: a soft blackbody with fcTg = 
68 eV and = 0.78A*(d/0.5)2 and a hard blackbody 
with fcTh = 320 eV and = 2.5 x 10-5^(^/0.5)2, 
where d refers to the distance in kpc. Here, A^ and A\^ 
indicate the observed emission erea of soft and hard 
blackbody components, respectively; A^ expresses the 
area of the whole neutron star surface. 

Mignani, Caraveo, and Bignami (1997) claimed 
that the blue/UV flux is consistent with an extrapo- 
lation of the soft, thermal part of the ROSAT X-ray 
spectrum. We thus compute the infrared flux from the 
Rayleigh- Jeans side of the soft X-ray spectrum with 
fcTs = 68 eV and A^ = 0.78A*(d/0.5)2. This IR pho- 
ton field is, in fact, so weak that the IC drag acting on 
a particle is negligibly small compared with the cur- 
vature drag. Therefore, the Vlasov equations can be 
solved without the IR field (i.e., the pair-production 
rate can be computed only from the GeV-keV coUi- 



3.2. Dependence on Injected Currents 

Let us now substitute the X-ray field into equa- 
tion (||) and solve the Vlasov equations by the method 
described in § ||. Once £'||(C)> '^±(0; and (7±*(0 are 
solved, we can compute the GeV and TeV spectra by 
the method described in § 4 in Paper VII. 

3.2.1. Equal Current Injection 

We first consider the case when ji = j2 . In figure ^ 
we present the acceleration field. 



ds 



UJp rUcC 



(28) 



for the four cases: ji — j2 — 0, 0.005, 0.05, and 2.0, 
which are indicated by the solid, dashed, dash-dotted, 
and dotted curves, respectively. In the third case 
for instance, the positronic and electronic currents 
equally fiow into the gap across the inner and outer 
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boundaries at the rate 0.05ri/27r per unit magnetic 
flux tube (i.e., 5% of the typical Goldreich- Julian 
value). For all the four cases, we set ai = 60°. Since 
the space charge is overfilled (or equivalently, —d^/ds 
distribution forms a brim as described in Hirotani & 
Okamoto 1998, see § 4.2 for details) if jgap exceeds 
typically several percent, we adopt jgap = 0.01 as 
the representative value. The abscissa designates the 
distance along the last-open field line and covers the 
range from the neutron star surface (s — 0) to the 
position where the disance equals s — O.SShjlc = 
3.29 X 10^ m. 

It follows from figure ^ that the gap is located 
around the null surface (where Bz vanishes), which 
is indicated by the downarrow. For a\ ~ 60°, 
n = 31.9 rad s"\ the nuU surface is located at s = 
1.59 X 10® m from the star surface along the last- 
open field line. The conclusion that the gap is located 
around the null surface for ji = j2, is consistent with 
what was obtained analytically in § 2.4 in Paper VII. 
We can also understand from figure |l| that W, and 
hence —d^/ds decreases with increasing ji = j2- This 
is because not only the produced particles in the gap 
(i-e., jgap) but also the injected ones (i. e., ji-\-j2) con- 
tribute to the 7-ray emission. Because of the efficient 
7-ray emission due to the injected particles, a sta- 
tionary pair-production cascade is maintained with a 
small pair-production optical depth W/Xp, and hence 
a small W. It is noteworthy that W does not lin- 
early decrease with jgap/jtot, by virtue of the 'nega- 
tive feedback effect' due to the factor; that is, 
the decrease in due to the reduced W partially 
cancel the original decrease in W. It may be note- 
worthy that iV^ cx WT cx W^^"^ holds (for details, see 
eqs. [8], [11], and [31] in Paper V, along with fig. 6 in 
Hirotani & Okamoto 1998). 

We present the 7-ray spectrum in figure |^; the 
curves correspond to the same cases as those in fig- 
ure ^ In each case, outwardly propagating 7-ray 
flux is depicted, because it is greater than the in- 
wardly propagating one. The observed fluxes are in- 
dicated by open circles and squares. It follows that 
the spectrum becomes soft as the injected current 
densities, ji = j2, increase. This can be easily un- 
derstood because —d^/ds decreases with increasing 
ji = j2 (fig. |l|). Even though solutions exist for su- 
per Goldreich- Julian currents (e.g., the dotted line in 
fig. H), the predicted spectra become too soft to match 
the observations. On these grounds, we can conclude 
that a sub- Goldreich- Julian current density is prefer- 



able for us to account for the observed spectrum. 

3.2.2. Current Injection across Outer Boundary 

We next consider the case when the gap is shifted 
towards the star surface by virture of the particle in- 
jection across the outer boundary. In figure we 
present —d'^/ds for vanishing ji] the solid, dashed, 
and dotted curves represent the cases of (ji,j2) = 
(0,0), (0,0.25), and (0,0.391), respectively. For all the 
three cases, we set a-i — 60° and jgap — 0.01. The 
abscissa covers the range from s = to s = 0.35n7Lc- 
There is no solution above j2 > 0.391, because the 
space charge is overfilled at the inner boundary above 
this critical value. 

It follows from the figure that the gap position 
shifts inwards as j2 increases. For example, the gap 
center position shifts from 0.225-n7Lc for j2 ~ (solid 
curve) to 0.049roLc and 0.0084ti7Lc for j2 = 0.25 
(dashed) and j2 — 0.391 (dotted), respectively, from 
the star surface along the last-open field line. In ad- 
dition, W decreases significantly as the gap shifts in- 
wards. This is because both the factors jgap/jtot and 
Ap decreases in equation (^7|). For details, see § 6.1 
in Paper VII. 

Let us now consider the Lorentz factor. As the 
maximum value in the gap, we obtain 3.0 x lO"^ for 
j2 — (solid curve), whereas 8.7 x 10® for j2 = 0.391 
(dotted one). That is, in the latter case, P reduces 
significantly even though the maximum of —d^/ds is 
greater than the former. This is because the parti- 
cles' motion becomes unsaturated due to the reduced 
W . Therefore, the mono-energetic approximation is 
no longer valid for the dotted curve in figure ||. 

The emitted spectra are presented in figure Be- 
cause the Lorentz factor decreases as j2 increases, the 
spectrum becomes soft as the gap shifts inwards. For 
the dotted curve, the spectrum would be further soft- 
ened from figure ^ if we considered the unsaturated 
effect of particle motion. We can therefore conclude 
that the gap should not be located well inside of the 
null surface for sa 60°, so that the predicted GeV 
spectrum may not be significantly inconsistent with 
observations. 

3.2.3. Current Injection across Inner Boundary 

Let us consider the case when the gap shifts to- 
wards the light cylinder, as a result of the particle 
injection across the inner boundary. In figure ^, we 
present —d'^/ds for vanishing j2] the solid, dashed. 
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B1 055-52 

inclination= 60 deg 




distance along field line [m] 



Fig. 1. — Distribution of the acceleration field, 
-d^/ds, for PSR B1055-52 for a-, = 60°. The 
solid, dashed, dash-dotted, and dotted curves 
correspond to the cases of (jgap,ji,j2)=(0.01,0,0), 
(0.01,0.005,0.005), (0.01,0.05,0.05), and (0.01,2,2), re- 
spectively. The downarrow indicates the position of 
the null surface, where vanishes. 



B1 055-52 

inclination= 60 deg 




distance along field line [m] 



Fig. 3.— Distribution of -d^/ds for PSR B1055-52 
for a-i = 60°. The solid, dashed, and dotted curves 
correspond to the cases of (jgap,ji,j2)=(0.01,0,0), 
(0.01,0,0.25), and (0.01,0,0.391), respectively. 
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Fig. 2. — Expected pulsed 7-ray spectra from 
PSR B1055-52 for a-, = 60° and jgap = 0.01. The 
four curves correspond to the same cases as presented 
in fig. 1. In each case, the larger flux for either the 
inwardly or the outwardly propagating 7-rays is de- 
picted. 



Fig. 4. — Expected pulsed 7-ray spectra from 
PSR B1055-52 for a-, = 60°, jgap = 0.01, and ji = 0. 
The four curves correspond to the same cases as pre- 
sented in fig. 3. In each case, the larger flux for either 
the inwardly or the outwardly propagating 7-rays is 
depicted. If the curve is thick (or thin) , it means that 
the outwardly (or inwardly) propagating 7-ray flux is 
greater than the other. 
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dash-dotted, and dotted curves represent the cases 
of (jij2)= (0,0), (0.25,0), (0.5,0), and (0.585,0), re- 
spectively. For aU the four cases, we set ai = 60° 
and jgap = 0.01. The abscissa covers the range from 
s = to s = l.Swi^c- There is no solution above 
ji > 0.585, because the space charge is overfilled at 
the outer boundary above this critical value. 

It follows from the figure that the gap shifts out- 
wards as ji increases. For example, the gap center 
is located at s = Scnt = {si + S2)/2 = 0.225n7LC 
for ji — (solid curve), while it shifts to the posi- 
tion s = 0.526n7Lc and s = 0.854n7Lc for ji = 0.25 
(dashed) and ji — 0.5 (dash-dotted), respectively. 
In addition, W increases as the gap shifts outwards. 
This is because the diluted X-ray field at the outer 
part of the gap increases Ap (see eq. Q). The in- 
crease of W with increasing ji, however, does not 
mean that —d'i'/ds increases with ji. This is because 
the small |pgj| oc \Bz\ at larger distance from 

the star results in a small \d{—d'^/ds)/ds\ in the Pois- 
son equation. 

The emitted spectra are presented in figure ^. 
The curvature spectrum becomes the hardest for 
(ji J2)=(0,0) of the four cases, because —d^/ds be- 
comes large by virtue of the strong magnetic field at 
relatively small distance from the star. It follows from 
the solid and dashed curves in figures ^ and || that 
the spectrum does not soften very rapidly when ji 
increases (fig. |) compared with the case when j2 in- 
creases (fig. |). This is because W increases with 
increasing ji (fig. ||), while it significantly decreases 
with increasing j2 (fig. ^). 

In passing, it is worth examining the relative mag- 
nitudes of the inward and outward fluxes. We con- 
sider the two cases ji = (solid) and 0.25 (dashed) 
for cki = 60°, jgap = 0.01, and j2 = in figure |^. 
The thick (or thin) lines denote the vF^, fluxes of out- 
wardly (or inwardly) propagating 7-rays at the outer 
(or inner) boundary. Because of the factor in the 
soft photon flux, the pair production is active at the 
inner part of the gap. Thus, the positron (or elec- 
tron) density rapidly increases (or decreases) with in- 
creasing ^ at the inner part, and nearly saturates at 
itot (or j2) in the outer part. See flgure ^ for the 
n+ (thick) and n-_ (thin) distribution for ai = 60°, 
igap = 0.01, ji ~ j2 = 0. Because of this asym- 
metric distribution of n_|_ and n_ along ^, outwardly 
propagating 7-ray flux dominates the inwardly prop- 
agating one. It is noteworthy that the extended gap 
for this middle-aged pulsar results in this large asym- 
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distance along field line [m] 

Fig. 5.— Distribution of -d^/ds for PSR B1055- 
52 for Qfi = 60°. The solid, dashed, dash- 
dotted, and dotted curves correspond to the cases of 
(jgap Jij2)=(0.01,0,0), (0.01,0.25,0), (0.01,0.5,0), and 
(0.01,0.585,0), respectively. 
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Fig. 6. — Expected pulsed 7-ray spectra from 
PSR B1055-52 for a; = 60°, jgap = 0.01, and j2 = 0. 
The solid, dashed, dash-dotted, and dotted curves 
correspond to the cases of ji — 0, 0.25, 0.5, and 
0.5847, respectively. Since the outward 7-ray flux 
dominates the inward one in each case, the former 
is depicted. 
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metry. For a small gap (i.e., W ^lc) obtained for 
young pulsars (Paper VII, VIII) and super-massive 
black holes (Hirotani & Okamoto 1998), the flux ra- 
tio between the outwardly and inwardly propagating 
7-rays are slightly less than unity, because and 
n_ distribute almost symmetrically with respect to 
the gap center. 

It also follows from figure ^ that the outwad 7-ray 
flux dominates the inward one when ji increases. This 
is because the injected positron at the inner boundary 
also contribute for the 7-ray emission. 

Let us summarize the main points that have been 
made in §§ p.2.l| - |3^2!3 . The curvature spectrum be- 
comes the hardest when ji — j2 — 0- In this case, 
the gap is located close to the null surface, where Bz 
vanishes. 

3.3. Dependence on Inclination 

It has been revealed that the 7-ray spectrum be- 
comes hardest when the injected current is small (say, 
when ji ~ j2 ~ 0) and that the EGRET spectrum 
cannot be explained for ai — 60° (as the solid line 
in figs. ^ or ^ indicates). Since the 7-ray energies 
are predicted to increase with increasing ai in Pa- 
per V, we examine in this section whether the EGRET 
pulsed flux around 6 GeV can be explained if consider 
a larger aj. In flgures ^ and ^ we present the accel- 
eration fields and the expected spectra for aj — 80°, 
75°, 60°, and 45° as the dash-dotted, dashed, sofid, 
and dotted lines. The current densities are chosen as 
jgap = 0.01, ji = j2 = 0. In each case, the outward 7- 
ray flux is depicted, because it dominates the inward 
ones. 



It follows from figure 10 that a large inclination 
angle (75° or greater) is preferable to explain the 7- 
ray spectrum around 6 GeV. The reasons are fourfold 
(see also § 6.2 in Paper VII): 

• The distance of the null surface from the star de- 
creases with increasing ai. 

• The magnetic field strength, and hence the Goldreich- 
Juhan charge density, QBz/{2Trc), in the gap increases 
as the distance from the star decreases. 

• It follows from the Poisson equation that the deriva- 
tive of —d'^/ds (i.e., —d^'^/ds^) increases with in- 
creasing Bz (i.e., with decreasing distance from the 
star). 

• By virtue of this increasing derivative, the accel- 
eration field, —d^/ds, at the gap center increases. 
We may notice here that W does not decrease very 
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Fig. 7. — Expected pulsed 7-ray spectra from 
PSR B1055-52 for a-, = 60°, jgap = 0.01, and = 0. 
The solid and dashed lines correspond to ji ~ and 
ji = 0.25, while the thick and thin ones to the out- 
wardly and inwardly propagating flux at ^ = ^2 and 
fi, respectively. 



Thick: positron flux Thin: electron flux 
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s [meter] 

Fig. 8. — Positronic (thick) and electronic (thin) den- 
sity distribution in the outer gap of PSR B1055-52 for 
tti = 60°, jgap = 0.01, and ji = j2 = 0. 
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incl.= 80,75,60,45° 




distance along field line [m] 



Fig. 9.— Distribution of -d^/ds for PSR B1055-52 
for ai = 80° (dash-dotted), 75° (dashed), 60° (sohd), 
and 45° (dotted). The current densities are fixed as 
jgap = 0.01 and ji = j2 = 0; therefore, the peak of 
each curve corresponds to the null surface for each 
inchnation. 
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Fig. 10. — Expected pulsed 7-ray spectra from 
PSR B1055-52 for 45° (dotted), a; = 60° (solid), 75° 
(dashed), and 80° (dash-dotted). The current den- 
sities are fixed as jgap = 0.01 and ji = j2 ~ 0. In 
all the four cases, outwardly propagating fluxes are 
greater than (or dominate) the inwardly propagating 
ones. 



rapidly with increasing ai, because of the 'negative' 
feedback effect due to the N-y factor in equation (p7|). 
As a result of this increased —d^/ds, the spectrum 
becomes hard for a large a^. 

It is worth noting that the Tev flux is always neg- 
ligibly small compared with the observational upper 
limit. This is because the infrared flux, which is de- 
duced from the Rayleigh-Jeans side of the soft surface 
blackbody spectrum, declines sharply as at small 
frequency, hv <C fcTg. 

3.4. Dependence on Distance 

We have assumed that the distance is 0.5 kpc so 
that the observed area of the soft blackbody emission. 
As = 0.78A,(d/0.5kpc)^, may be less than A,. In 
this section, we consider the case when the distance is 
1.5 kpc, as indicated by Taylor, Manchester and Lyne 
(1993), and compare the results with d = 0.5 kpc case. 

In figure |ll|, we present the resultant spectra for 
ai = 80°, jgap = 0.01, and ji = ja = 0. The solid (or 
dashed) line represents the 7-ray spectrum when the 
distance is 0.5 kpc (or 1.5 kpc). The normalization 
is adjusted by the fluxes below 0.5 GeV; namely, the 
cross sectional area is D\ = 2.3^VF^ (or 7.7^PF^) is 
assumed for d = 0.5 kpc (or 1.5 kpc). It is clear from 
figure |ll| that the GeV spectrum becomes hard if we 
assume a smaller distance, d. 

The conclusion that a small distance to this pul- 
sar is preferable, is derived, because the less dense 
X-ray field for a smaller d results in a larger Ap in 
equation (|2^), and hence in a larger W and the ac- 
celeration field, —d'^/ds. However, the solution does 
not vary so much even though the X-ray density in- 
creases (1.5kpc/0.5kpc)^ = 9 times. This is due to the 
'negative feedback effect' caused by the factor in 
equation (p7|); that is, the increase of W is substan- 
tially canceled by the increase of N^^ which is propor- 
tional to the product of the curvature-radiation rate 
of a particle per unit length and W . In other words, 
the solution exists for a wide range of parameters and 
does not change very much if the variation of the X- 
ray field density is within one order of magnitude. 

4. Discussion 

In summary, we have developed a one-dimensional 
model for an outer-gap accelerator in the magne- 
tosphere of a rotation-powered pulsar. Solving the 
Vlasov equations that describe a stationary pair- 
production cascade, we revealed that the accelerator 
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shifts towards the star surface (or the hght cyhnder) 
as the particle injection rate across the outer (or in- 
ner) boundary approaches the Goldreich-Juhan value. 
Applying this theory to a middle-aged pulsar, B1055- 
52, we find that stationary solutions exist for a wide 
range of parameters such as the injected currents, 
magnetic inclination, and the distance to the pulsar. 
Comparing the expected spectrum with ERGET ob- 
servations, we conclude that a sub Goldreich- Julian 
current density, a large magnetic inclination, and a 
small distance (500 pc, say) are preferable for this 
pulsar. 

4.1. Created current as a Free Parameter 

In this subsection, let us discuss the physical rea- 
son why jgap should be treated as a free parame- 
ter, rather than a quantity that is solved from the 
Vlasov equations. We first consider the vacuum limit 
in which the true charge density {B /Bcn_t){n'+ — n ) in 
the Poisson equation (||) is negligible compared with 
the Goldreich- Julian charge density, Bz/Bcnt- In this 
limit, it follows from equations (|l^) and (^4|) that if 
n± and g^. form a solution of the Vlasov equations 
then /an± and /a5± do, as long as the arbitrary fac- 
tor /a satisfies /a[n+ — n-] <C Bz/B. It is jgap that 
parameterizes the arbitrary factor, /a. 

Next, we consider a non- vacuum gap in which the 
screening effect due to [B / Bcnt){n+ — n_) becomes 
important in the Poisson equation. In this case, the 
solution deviates from the vacuum one as a continu- 
ous function of jgap (see figs. 9, 10 in Paper I, where 
Jo in that paper corresponds to jgap in the present pa- 
per). The same can be said when the inverse Comp- 
ton scatterings dominate (figs. 8, 9, 10 in Paper II; see 
also eq. [50] , which describes how the voltage drop de- 
pends on jo, i.e., jgap). The nature of jgap as a free 
parameter is unchanged if particles are injected from 
the boundaries. For example, equation (^) shows 
that the pair production physics determines the gap 
width as a function of jgap/jtot, even if ji and j2 are 
fixed. 

In the present paper, we are primarily interested in 
the outer-magnetospheric accelerator. Therefore, to 
emphasize that a gap itself has a freedom, we regard 
jgap as a free parameter, in addition to ji and ^2, and 
calculate jtot from jtot = jgap + ii + 32- However, 
we may alternatively consider that the total current 
jtot is determined globally (e.g., in the wind region). 
In this case, we may have to regard jtot, ji, and ^2 
as free parameters in an outer gap and calculate jgap 



from jgap = itot-ji-j2- In another word, jgap may be 
constrained by a global condition, if we consider the 
polar-cap, outer-magnetospheric, and wind regions as 
a whole. 

4.2. Maximum Current Density 

We next discuss why there is an upper limit for 
the current flowing in the gap. It follows from the 
Poisson equation (^) that dE\\/d£, changes its sign 
where a space charge is overfilled in the sense that 
|n4. — n_| exceeds \Bz\IB, or equivalently, in the sense 
that jgap/odd(a;) exceeds (ji - j2)feven{x) + Bz/B in 
equation (30) in Paper VII, where x designates the 
position in the gap (see § 2 of Paper VII for details) . 

If the 'overfilling' occurs at a certain position, 
dE^^/d^ changes its sign there and will diverges out- 
side of the gap (fig. 2 in Hirotani & Okamoto 1998). 
Under this circumstance, we obtain two accelerators: 
the original one and an additional one outside of the 
point where is minimum. (Usually, such an over- 
filling occurs close to one of the two boundaries, where 
|n_|_ — 71-1 niaximized.) Since diverges in the lat- 
ter accelerator, a large fiux of 7-rays are expected to 
illuminate the original gap. The resultant, excessive 
7-ray flux leads to an additional pair production in 
the gap, which screens the original This screen- 
ing may continue until closes at the boundaries 
and hence until the additional accelerator vanishes. 
We conjecture that such a solution having diverging 
(and hence belonging to a higher energy state) is 
dynamically unstable. We thus consider that there 
is an upper limit for Jgap 

above which stationary so- 
lutions cease to exist due to the overfilling of space 
charges. 

4.3. Pulse Phases at Different Frequencies 

We next consider the expected pulse-phase rela- 
tionship between the incoherent GeV emission from 
the outer gap and the coherent radio emission from 
the polar cap (e.g., Michel 1991; Zhang & Harding 
2000). If the neutron star is nearly orthogonally ro- 
tating (ofi > 75°), as indicated in § |3^ (or fig. and 
if the outer gap is located near to the null surface, 
as indicated in the last paragraph n § 3.2, both the 



polar-cap and outer-gap emissions are beamed nearly 
in the same direction (fig. ^2|). It is noteworthy that 
the intersection of the null surface and the last-open 
field line is located only at 0.0278tZ7LC (or equiva- 
lently, 8.71 X 10~*s in light crossing time) from the 
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star for aj = 80°. It is therefore possible that the for- 
mation of an outer gap near the null surface for a large 
inclination angle explains the observational fact that 
the GeV peak in the light curve occurs at the close 
phase with one of its two radio pulses (Thompson et 
al. 1999, and references therein). 

4.4. Spectrum at Different Positions 

In this subsection, we present the 7-ray spectrum 
as a function of the position in the gap. As the best- 
fit case, we choose = 80°, ji = j2 — 0, and 
igap = 0.01. To integrate the 7-ray flux in certain 
regions in the gap, we divide the gap into the follow- 
ing four regions: 

1.006 + 0.006 < < 0.75a + 0.256, 
0.756 + 0.256 < <e < 0.506 + 0.506, 
0.506 + 0.506 < ^ < 0.256 + 0.756, 
0.256 + 0.756 < < 0.006 + 1-006- 
^, we present the outward 7-ray flux emit- 
ted from each region: The solid, dashed, dash-dotted, 
and dotted lines correspond to the fluxes from re- 
gions 4, 3, 2, and 1, respectively. It follows from the 
figure that the outward 7-ray fluxes from regions 2, 
3, and 4 do not differ very much. However, from re- 
gion 1, the flux is much smaller than those from other 
regions. The figure also indicates that the spectrum 
becomes hardest from region 3, because the electric 
field is strong there. The 7-rays emitted from region 1 
becomes soft, mainly because the curvature radius is 
large at the inner region. 

4.5. Magnetospheric Current Distribution 

As demonstrated in § 6.3 in Paper VII, the spin- 
down luminosity of a rotation-powered pulsar be- 
comes (see § 6.3 in Paper VII) 



region 1 
region 2 
region 3 
region 4 

In figure IT 



Er, 



' jtotfi 



2 

active 



(29) 



where fim is the neutron star's magnetic dipole mo- 
ment; /active represents the ratio of magnetic fluxes 
along which currents are flowing and those that are 
open. ForPSRB1055-52,17Vm/c^ = lO^^'-^^'ergs s'^. 
If all the open field hues are active (i.e., /active = 1), 
Jtot ~ 1 is required so that the observed spin-down 
luminosity 10^'^'^^ergs s~^ may be realized. If a small 
fraction of open field lines are active (i.e., /active 1), 
a super-Goldreich- Julian current (i.e., jtot 3> 1) is re- 
quired. 

However, our outer-gap model requires (and, in 
fact, previous outer-gap models assume) jtot 1- 
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Fig. 11. — Expected pulsed 7-ray spectra from 
PSR B1055-52 for a; = 80°, jgap = 0.01, and ji = 
j2 = 0. The solid (or dashed) line represents the 7-ray 
spectrum when the distance is 0.5 kpc (or 1.5 kpc). 



polar-cap emission 
^ outer-gap emission 



outer gap located close 
to the null surface 




Fig. 12. — Side view of an orthogonal rotator. Radio 
emissions are supposed to be emitted from the polar 
cap, while the 7-rays are from the outer gap. 
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Thus the deficient current must flow along the open 
field lines that are not threading the accelerator. 

4.6. Comparison with Previous Works 

Let us finally compare the present work with the 
numerical simulation performed by Smith et al. (2001) , 
who examined equilibrium charge distributions in the 
magnetosphere of an aligned rotator. In the first 
place, they demonstrated that the Goldreich- Julian 
charge distribution is unstable and collapses to form 
a polar dome containing plasma of one charge and an 
equatorial belt containing plasma of the other sign: 
E • B = inside both of them. These are separated 
by a vacuum gap in which E • B 7^ holds. To apply 
the theory developed in the present paper to the case 
when the Goldreich-Julian charge distribution breaks 
down, we must modify the Poisson equation (eqs. Q 
and [eq:basic-2]) so that the electric field caused by 
the charges in the dome and the belt may be taken 
into account. 

In the second place, Smith et al. (2001) demon- 
strated that the creation of electron-positron pairs in 
the gap between the dome and the belt reduces the 
value of E • B in the gap so that it turns off at last. 
This is because the charges are attracted to regions of 
the same sign, following the closed field lines, and in- 
creasingly filling the magnetosphere. Instead, as con- 
sidered in the present paper, if pairs are created in 
the open field line region, we may expect global flows 
of charged particles and hence a stationary pair pro- 
duction cascade in the magnetosphere of a spinning 
neutron star. If we extend the work by Smith et al 
(2001) into oblique rotators and apply the present 
method to the gap, we may construct a stationary 
magnetospheric model with pair production cascade 
under the existence of global currents. There is room 
for further investigation on this issue. 

One of the authors (K. H.) wishes to express his 
gratitude to Drs. Y. Saito and A. K. Harding for valu- 
able advice. He also thanks the Astronomical Data 
Analysis Center of National Astronomical Observa- 
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Fig. 13. — Gamma-ray flux from PSR B1055-52 from 
different positions in the gap, for cci = 80°, jgap = 
0.01, and ji = j2 = 0. The solid, dashed, dash- 
dotted, and dotted lines represent the vF^ flux from 
regions 4, 3, 2, and 1, respectively (see text). 
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